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ABSTRACT 

Heavy nuclei such as nickel-56 are synthesized in a wide range of core-collapse supemovae (CCSN), in- 
cluding energetic supemovae associated with gamma-ray bursts (GRBs). Recent studies suggest that jet-like 
outflows are a common feature of CCSN. These outflows may entrain synthesized nuclei at launch or during 
propagation, and provide interesting multi-messenger signals including heavy ultra-high-energy cosmic rays. 
Here, we investigate the destruction processes of nuclei during crossing from the stellar material into the jet 
material via a cocoon, and during propagation after being successfully loaded into the jet. We find that nuclei 
can survive for a range of jet parameters because collisional cooling is faster than spallation. While canonical 
high-luminosity GRB jets may contain nuclei, magnetic-dominated models or low-luminosity jets with small 
bulk Lorentz factors are more favorable for having a significant heavy nuclei component. 
Subject headings: gamma-ray burst: general — nuclear reactions, nucleosynthesis, abundances supemovae: 
general 



1. INTRODUCTION 

Recent studies have led to a canonical picture where 
long gamma-ray bursts (GRBs) are rare types of core- 
collapse supemovae (CCSN) that are acco mpanied by the 
launch of energetic relativ istic jets (see, e.g.. lMeszarosll20()6c 
IWooslev & BloomI 12006'. for reviews). While the GRB is 
a rare phenomenon, a significantly larger fraction of CCSN 
could produce jets that do not successfully produce GRBs 
either because of energetic or collimation reasons. Many 
such jets may even be choked in their progenitor envelopes 
dMeszaros & WaxmanI l200lb . Observationally, spectropo- 
larimetry shows that the degree of asymmetry in CCSN in- 
creases with time and hence to greater depths of the CCSN, 
suggesting an association with the central engine and possi- 
ble la unch of bipolar jets (IWang et al.ll200ll: IChomock et al.l 
I2OIOI) . SN 1987A, the closest CCSN in modem times, shows 
a globally asymmetric expanding d ebris with an axi s that 
roughly aligns with that of its rings (IWang et al.ll2002ft . Cas 
A, one of the well-known CCSN remnants in our Galaxy, may 
have been accompanied by an iron-rich jet component which 
may help t o explain the observed apparent overturn of the Fe- 
rich ejecta dWheeler et al.ll2008t lOeLanev et al.ll2010l) . 

It is well known that nuclei can be synthesized not only 
in stellar nucleosynthesis but also during the CCSN, where 
pre-existing nuclei in the star are further fused into heav- 
ier nuclei via explosive nucleosynthesis. For CCSN as- 
sociated with GRBs, such as SN 2003dh (GRB 030329) 



and SN 1998bw (GRB 980425), lar-ge kinetic energies of 
^ 10^^ erg are suggeste d and the inferred synthesized ^^Ni 
masse s are ^ 0.5 M p-, flwam oto et al.|[l998t IWooslev et alj 
199f^; IWooslev & He ger 2003). The large explosion ener- 
gies may be due to a baryon-rich jet component driving 
the CCSN, which would lead to different elemental yields 
due to additional nucleosy nthesis and/or efficient mixing and 
transport of stellar nuclei (Lemoine 2002; P met et all 120021 : 
iMaeda et alj|2002l iMaeda & NomotQ,2003.) . Alternatively, 
a large amount of ^^Ni can be produced by a disk wind, 
where free n ucleons ejected from th e disk combi ne to form 
heavy nuclei dMacFadven & Wooslev 1999; Pmet et al.ll2003l 



ISurman & McLaughlinluOOSb . The highly relativistic jet re- 
sponsible for the GRB should have a smaller baryon mass, 
but may still pick-up and/or entrain stellar nuclei, nuclei syn- 
thesized by a wider jet, or nuclei synthesized by a disk wind. 
The fate of nuclei in relativistic jets is also interest- 
ing in light of recent reports on the nuclear composition 
of ultra-high-energy cosmic rays (UHECRs). The origin 
of UHECRs remains one of thegreat mysteries in high- 
energy astrophysics (see, e.g., JHil las 2005; Bliim er et al.l 
l2009t iBeattv & Westerhoil2009t iKotera & Olintoll201ll for 



reviews). At the highest energies (above '^ 50 EeV), 
the extreme energies argue for extragalactic sources (al- 
though Galactic transients may give a contribution around 
^ 10 EeV, e.g., Calvez et all |2010|) . Source candidates 
fall into active galactic nuclei (AGNs), including radio- 



loud AGNs dBiemiann & Strittmatted [ l987t lTakahara|ll990[ 



lNormanetal .i l 1995b. radio-quiet AGNs (iPe'eret all (20091) 
young AGNs ('Takahara" 1 993 |Takami & Horiuchi"201 1), as 
well as powerful tr ansient fla res (iFarrar & Gruzinov 2009, 
Dermer et al.l l2009l): G RBs JWaxmanl 1 19951; IVietril 1 19951 



Mil^rom & Usovlll995b . including sub-luminous GRBs ac- 



companying relativistic ejecta in some form (Murase et aL 



gOOetlWang et al.l2007l;jMurase et al.l2008l;IChakraborti et al 



201 it iLiu & WangI 1201 lb : the formation of rapidly ro- 
tating strong l y mag n etized protorn a gnetar s dAr ons 2003; 
Murase etal.ll2009t iMetzger et all boiU Koteral 120111; 
Fanget al. 2012b; and galaxy cluster sh ocks dNorman et al.l 



1995l;lKang et alJll996l;[lnoue et alj2007l) . Non-astrophysical 



sources include top-down scenarios where the UHECR arise 
from decays of massiv e relics from the early universe (e.g., 
lBerezinskvetal.lll997b . 

The composition of UHECRs is observationally inferred 
by measuring composition-dependent quantities of the show- 
ers made as UHECRs enter Earth's atmosphere. The Pierre 
Auger Observatory (PAO) finds that the average shower depth 
at shower maximum, Xmax, and its rms variation rms(X,„ax), 
suggest the composition becoming increasingl y dominated by 



heavy nuclei above the 'a nkle' of lO^'*-^ eV (lAbraham et al.l 
|2010[ lAbreuetal.ll2011ab . Note that the PAO results have 
not be en verified by the H igh Resolution Fly's Eye experi- 
ment flAbbasi et al.l 2010ab and th e Telescope Array experi- 
ment (e.g.. iTsunesada et al.ll201 1|) . although the latter exper- 
iments observe a different hemisphere and their statistics are 
lower. In addition, these experiments do not confirm the di- 
rectional anisotropies reported by PAO (lAbraham et al.ll2007l 



lAbreu et al.ll2oTot lAbbasi et allbOOSl bOlObV Also, the PAO 
indicators have been sho wn to be inconsistent with each other 
dShaham & Piranll2012b and that they may not reflect a com- 
position change given the current uncertainty in shower in- 
teraction physics ( Wilk & Wlodarczvk 2011). More observa- 
tional and theoretical studies are required to settle these and 
other remaining issues. 

If UHECRs are indeed composed of heavy nuclei, this pro- 
vides important implications for UHECR sources. In the 
AGN and galaxy cluster shock origin scenarios, the matter 
being accelerated originates from the intergalactic medium, 
so that the mass fraction of Fe and heavier nuclei is small: 
^ 10""^ for solar metallicity. Although the dominance of nu- 
clei in UHECRs may be explained by a rigidity-dependent ac- 
celeration mechanism, where the maximum Fe energy would 
be Zpc = 26 times higher than those of protons, this requires 
the maximum accelerated energies in all contributing sources 
to be somewhat fine-tuned such that fi'p^max ~ 10^^ eV. 
Also, such a scenario predicts a larger pr oton-to-nuclei ratio at 
low energies than is actual ly observed (ILemoine & WaxmanI 
120091; lAbreu et ail l2011bb . An alternate scenario is that 
the environment supplying the injected particles is enriched 
in nuclei. The abundance of nuclei inferred from X^ax 
and rms(Xmax) require high nuclei abunda nces (see, e.g., 
lAnchordoqui et aU 120071; iTavlor et akl l201ll and references 
therein). Massive-star origins, including GRBs, CCSNe with 
relativistic ejecta, and magnetars would be attractive sources 
in this sense, because as described above the environment 



contains large fractions of intermediate or heavy nuclei. Ad- 
ditionally, it has been shown that once loaded, nuclei may 
be accelerated to ultra-high energies and successfully sur- 
vive in the dissipation regions of jets, including both classical 
high-luminosity GRBs (HL GRBs) a nd sub-luminous GRB 
dMurase et al.ll2008l; IWang et al.ll200 8V 

In this paper, we investigate the origins and survival of nu- 
clei as the jet is still inside the star More specifically, we 
first discuss locations where nuclei may enter the jet medium 
and investigate whether nuclei can survive in each of these 
locations. Second, we investigate whether nuclei that have 
successfully made their way into the jet survive during the 
jet propagation through the progenitor star These issues are 
different from previous works whose main focus were on the 
survival of UHECR nuclei in the emission regions (where the 
GRB occurs, typically outside the star). Note that UHECR 
acceleration is not expected inside the star even though our 
study would be useful for UHECR sources, and we mainly 
consider the survival of low-energy nuclei. 

The paper is organized as follows. In Section|2] we discuss 
sources of nuclei and discuss conditions for nuclei survival in 
each of them. In Section |3j we consider processes in the jet 
that may lead to nuclei disintegration and discuss conditions 
under which nuclei survive. Finally, we finish with discus- 
sions and conclusions in Section |4| We express quantities as 
Qx = Q/IO^ in cgs units. 

2. NUCLEI ORIGINS AND SURVIVAL 

We first discuss three potential sources of nuclei in the con- 
text of GRBs. First, loading at the jet base, where a small 
amount of baryons may be loaded into the jet (the jet is as- 
sumed to be radiation or magnetically dominated). Depending 
on the jet parameters, nuclei survive or are disintegrated into 
free nucleons. Second, nucleosynthesis in the outflow itself, 
which occurs if nuclei are disintegrated into free nucleons at 
initial loading. Finally, entrainment, which is the loading of 
surrounding nuclei into the jet during jet propagation. 

2.1. Nucleus Loading at the Jet Base 

The progenitors of GRBs are massive stars with Fe cores of 
mass approximately 1-2^/© that extend to radii a few x 10* 
cm. It is thought that a few seconds after the onset of col- 
lapse, a compact object (either a neutron star or a black hole) 
forms at the center, surrounded by an equatorial disk sup- 
ported by centrifugal forces. By some currently unconfirmed 
mechanism, an energy reservoir is tapped and energy is re- 
leased near the collapsing core. In the "collapsar" model, the 
energy derives from the gravitational energy of rap id accre- 
tion, and neutrinos play the role of en ergy transport (IWooslevI 
119931; iMacFadven & WooslevI [T999h . In this case, the en- 
ergy is deposited dominantly as radiation energy. Alterna- 
tively, the jet may be dominated by a magnetic component, as 
for example in scenarios where the jet magnetic field threads 
the black hole event horizon and the power derives from the 
black hole spin energy ( Blandford & Znajek 19 77), where the 
j et magnetic field threads the surface of the accretion disk 
dBlandford & Pavnd [19821 iProga et al.1 l2003b . or where the 
power derives from the spin-down energy loss of a central 
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Fig. 1. — Optical depths for photodisintegration (solid) and spallation 
(dashed) as functions of jet injection radius, for various jet radiation 
luminosities. For a pure radiation fireball, Lrad.o corresponds to Lj^^ and 
the GRB luminosity by Equation (34). In the case of a magnetic-dominated 
outflow, Lrad.o ^^Y t>e significantly smaller for a given observed GRB 
luminosity, and nuclei loading would be easier Bottom panel shows the 
dominant pre-supemova stellar abundances taken from the rotating 2OAf0 
pre-superaova progenitor model E20 of Heger et al. (2000). The abundances 
are illustrative: depending on the timing of jet launch, nuclei at small radii 
will be disintegrated by a supernova shock or synthesized by explosive 
nucleosynthesis before the jet propagates through (see the text). 



magnetic protoneutron star (lUsovll 1992b . 

We consider a jet with baryonic matter injected at rate Mq, 
radiation energy at rate irad.o ^ Mqc^, and magnetic energy 
at rate L^ag,o (we will use coUimation-corrected values, not 
the isotropic equivalent) at a radius vq. The jet subsequently 
adiabatically expands as its radiation and/or magnetic ener- 
gies are gradually converted to bulk kinetic energy. The max- 
imum bulk Lorentz factor is set by 77 = Lkc/(^^c^), where 
Lkc = .^rad,o + imag.o IS the kinctic luminosity after the end 
of the jet bulk acceleration phasqj. We defined a "fireball" as 
a jet where radiation dominates, i.e., Liad.o ^ -^mag.o^ and a 
magnetic outflow as the opposite with imag.o ^ Lrad.o- 

The large values of 77 > 100 required in GRB (iRudermanI 
1 1975b necessitate a jet with small baryon loading. The quan- 
titative predictions about where and how baryon loading oc- 
curs are tentative at best, but any baryons must come from the 
surrounding material. The composition of the pre-supemova 
stellar core is dominantly nuclei (bottom panel of Figure [1}. 
At small radii this can be significantly altered by the col- 
lapse. For example, if a supernova shock is launched prior 
to the jet, nuclei out to a distance of ^ 10^-10^ cm would be 
disintegrated by the supernova shock. At later times, explo- 
sive nucleosynthe sis can alter the composition. For example, 
iFrver et al.l (l2006b demonstrate with one-dimensional simula- 

' We caution that in magnetic-dominated jets such definitions may not 
strictly hold because of the model-dependent conversion of magnetic energy 
to kinetic energy; see Section |2.3.1| 
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Fig. 2. — Constraints on jet parameters for Fe nuclei survival during load- 
ing at the jet base. Photodisintegration (solid) and spallation (dashed) are 
shown separately for several values of rg. Shaded regions result in nuclei 
destruction, while white areas allow nuclei survival (shades shown for our 
nominal tq = 10^ cm). For magnetic-dominated jets, the situation is more 
favorable and quantitatively depends on the jet energy composition (see the 
text). 

tions that up to ~ IMq of ^^Ni may be synthesized, depend- 
ing on whether the core collapse proceeds to a black hole di- 
rectly or via fall back, and also on the explosion energy. The 
detailed composition of the surroundings therefore depends 
on the GRB scenario and timing of jet launch relative to core 
collapse. In either case, it is expected that there are abundant 
nuclei especially for radii larger than ^10^ cm. 

Any nuclei that are loaded into the jet are destroyed if it 
collides with particles or photons with energies exceeding the 
nuclear binding energy ^ 10 MeV in the nuclei rest frame. 
We consider a jet of luminosity irad,o injected at a radius rg 
with initial Lorentz factor Fq = 1. Its temperature is related 
to the radiation energy density. 



aT^ 



-^rad,0 



(1) 



where a is the radiation constant and Sq ~ f^o^o^ is the jet 
cross section. Typically, the jet is collimated already at injec- 
tion due to rotation or magnetic geometries, and we adopt a 
solid angle Og =0.1 sr independent of tq. This results in 
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-■^■■J "0,-1 ^rad, 50^^0,7 J^Vie V . 



(2) 



Therefore, nuclei will be destroyed by the high energy tail of 
the photon spectrum. The optical depth for photodisintegra- 
tion is TA-i — n-fdA-yfo/^Q, where n^ is the photon density 
and (TA-y is the photodisintegration cross section which peaks 
at approximately 10"^^ cm'^ (e.g.. iMurase et al.ll2008b . We 
require that ta-/ < 1 for the survival of nuclei during loading. 
Similarly for spallation we require that Tgp ~ nQaspro/To < 
1, where no is the comoving ion density. 
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and A is the average mass number of the baryons in the jet. 
Also, asp = oqA^I'^ with ctq ~ 3 x lO^^e ^^ ^^^ A = 56 
is the mass number of the injected nuclei of interest (assumed 
Fe). 

In Figure [T] the photodisintegration and spallation optical 
depths are shown for three values of Lrad.o- We adopt a 
thermal photon spectrum of temperature Tq, and assume that 
the jet ions have a Maxwell-Boltzmann velocity distribution 
with temperature equivalent to the photon temperature. Pho- 
todisintegration is more important than spallation, because the 
larger photon density compensates the slightly smaller peak 
photodisintegration cross section. Note that in Figure [T] we 
adopt j4 = 1, which gives the largest spallation target density 
possible. 

As expected, for a GRB fireball of typical luminosities 
(L,ad,o = lO'^^-lO^iergs-i), the high Tq results in nuclei 
being photodisintegrated at loading. However, since photo- 
disintegration occurs with the exponential tail of the photon 
spectrum, nuclei survival is highly sensitive on Tq and hence 
GRB parameters. For example, nuclei survival is possible 
for To greater than 10^-10^ cm. Values beyond 10^ cm are 
somewhat large but plausible if the jet is powered by an accre- 
tion disk. However, in this case it may be difficult to achieve 
short variable time scales on the order of ms observed in some 
GRBs, unless they are imprinted by instabilities during jet 
propagation. Figure [1] also demonstrates that nuclei survival 
is more easily possible for sub-luminous GRBs. We quantify 
these constraints in Figure |2] 

It is easier to initially load nuclei in magnetic-dominated 
jets (£mag,o ^ iiad,o); bccausc £rad,o IS Smaller for a given 
observed GRB luminosity and Tq is correspondingly smaller 
Both photodisintegration and spallation become less effective. 

2.2. Nucleosynthesis by Outflows 

The fireball is composed of free nucleons if the photodisin- 
tegration optical depth at jet launch is larger than unity. As the 
jet expands and cools, nucleons combined to form nuclei. The 
freeze-out composition of HL GRB jets has been investigated 
by various authors, and strongly depends on the radiation en- 
ergy content of the jet. For a canonical GRB fireball of initial 
temperature ^ 1 MeV, the entropy per baryon is necessarily 
high, S > 10^ kb nucleon^^. Free nucleons combine to form 
a particles only once the deuterium bottleneck has been bro- 
ken, which for higher entropies occurs at lower jet densities 
where the further processing into carbon and higher nuclei are 
not fast enough compared to the expansion time scale. As a 
result, few elements heavi er than He are f ormed, similar to 
big bang nu cleosynthesis (iLemoind 120021; iPruet et akl 120021; 
lBeloborodov.2003.) . 

The situation is dramatically different for magnetically 
dominated jets, where most of the energy is stored in mag- 
netic energy and the jet entropy is necessarily lower, S ^ 10- 
300 fcfc nucleon"^. Under such conditions a recombination 
occurs at higher densities such that heavier elements can be 
formed efficiently via the triple-a process and subsequent a 
captures. Focusing on the magnetic jets of a mil lisecond pro- 
tomagnetar central engine, iMetzger et al.l (120111) showed that 
the jet composition may indeed be dominated by heavy nu- 



cle i. Also, depending on how neutron-rich the jet matter 
is (iMetzger et al.ll2()()8h . neutron-catalyzed a formation and 
neutron captures open the path for the rapid synthesis of Fe- 
peak and heavier n-capture nuclei (e.g., A > 90). Similar 
nucleosynthesis i s also realized in m ore baryon-rich jets with 
lower entropies (llnoue et al.ll2003l) . Although the baryonic 
jets cannot produce the classical GRB phenomenon, it can be 
relevant for sub-luminous GRBs and hypernovae where rj can 
be much lower than rj ^ 100-1000. 

When the central engine consists of a black hole and ac- 
cretion disk system, hot outflows from the disk are expected, 
where the disk wind with S ~ 10-100 kf, nucleon"^ can ac- 
count for the amount of ^^Ni observed in hypernovae asso- 
ciated with GRBs (e.g.. ISurman & McLaughlinll2005h . We 
discuss the entrainment of external nuclei next. 

2.3. Nucleus Entrainment During Propagation 

As the jet propagates through the star, it can pick up 
baryons from the stellar core, surrounding wider jet, or the 
disk wind environment. Since these environments can be 
nuclei-rich, we discuss how nuclei may survive during the en- 
trainment process. We will assume a proton-dominated jet 
(A — 1) for this purpose. For our canonical jet, we con- 
sider a classical high-luminosity GRB jet, fixing Fq = 1 
and 57o = 0.1 sr. Further, we adopt Lkc — 10^°ergs~^, 
rj = 100, and rg = 10^ cm, and show dependencies where 
applicable. For illustration purposes, we adopt the rotating 
20Mf7) pre-supernova progenitor model E20 of iHeger et al.l 
(I2OOOI) . whose density profile can be well approximated by 
p*(r) = 1.7 X 105(r/10^cm)-"gcm-3withn«3. 

2.3.1. Cocoon Properties 

After the onset of collapse, the lack of centrifugal force 
along the rotational axis leads to free mass infall and the for- 
mation of a funnel region where conditions are favora ble for 
the launch of a jet ( Wooslev 1993; MacFadven & W ooslevI 
Il999h . In magnetic-dominated jets, the outflow initially flows 
along open field lines. However, self-collim ation via mag- 
netic hoop stress fails in relativistic flows (e.g.. lSakuraill 19851; 
'Bucci antini et al.l I2OO6 ). and recent works show how inter- 
actions with the star may redirect the flow toward the poles 
and produce relativistic bipola r jets along the rotational axis 
(Bucciantini et al. '2007', 2008; iKomissarov & Barkovl I2OO7I; 
.Uzdenskv & MacFadven.2007i) . We model the jet cross sec- 
tion as 

^{r) = ^J^\ , (4) 

where for example ^ = 2 corresponds to a conical (or spheri- 
cal s ymmetric) jet, and £ = 1 corresponds to a funnel-like jet 
(e.g.,'Meszai-os & Rees 2001"). 

After injection, the GRB jet accelerates as its radiation or 
magnetic energy is converted to bulk kinetic energy. For a 
radiative fireball expanding adiabatically wi th cross section E, 
energy conserva tion yields (see also, e.g., iToma et al.,i2007l ; 
llokaetalJl201lh 



r,M 



ro(r/ro)«/2 



r < rsat 



(5) 



where rsat = V^^^^q is the saturation radius where the ter- 
minal Lorentz factor is reached. Often r < Vgat is called the 
jet or bulk acceleration phase and r > rsat is the jet coasting 
phase. 

The relativistic jet slows down abruptly in a narrow layer at 
the head of the jet where it comes into contact with the over- 
laying stellar material. Due to the large ram pressure experi- 
enced, a strong reverse shock forms and the jet head is decel- 
erated to a sub-relativistic velocity /?/,<!. Further upstream, 
a forward shock forms that shock heats the stellar material. 
The jet head velocity is set by balancing the ram pressures ap- 
plied on the forwa rd and reverse s hocks. We use the analytic 
approximations of lMatznerl (120031) . 



Ph = P, 



1 



l + Z-l/2' 



where L is a dimensionless parameter defined as 



L- 



^{r)p.c^' 



(6) 



(7) 



In the limit L -^ I, which corresponds to a relativistic reverse 
shock (as viewed in the jet frame) as is generally the case for 
jet propagation in the star. Equation (|6]l can be linearized, and 



Ph 



L''^ oc ikc^/V("-«)/V(«-^)/^ 



(8) 



We call the shocked stellar and shocked jet material collec- 
tively as the jet head. At first, the jet head pressure is smaller 
than the surrounding stellar pressure and the jet head remains 
pressure confined. Thejet head pressure is Ph — C//i/3, where 
Uh — ^Tj^npTTipC^ is the jet head energy density and Up is 
the proton density, i.e., the ion density 



n,,(r) 



-^ko 



T,{r)rjT j ArUp 



(9) 



with A = 1. Hence Ph falls as ex r^^' ^ during the bulk accel- 
eration phase and ex r~^ during the coasting phase, while the 
stellar pressure falls as P, ex pj^^'^ ex r^^n/s Therefore, the 
jet head pressure eventually overtakes the stellar pressure and 
starts to overflow. We call the radius at which the pressures 
become equal as rc\ cocoon formation occurs for r > Tc- 
For typical parameters, Vc is in the range 10^-10^ cm and 
cocoon formation occur s during the bulk acceleration phase 
dMeszaros & Reesll200l]) . 

For magnetic-dominated outflows, the conversion of mag- 
netic energy to kinetic energy is prolonged since only part 
of the magnetic luminosity gets converted directly into 
kinetic luminosity, the other part being converted to ki- 
netic energ y in a two-stage proces s through thermal en- 
ergy (e.g., iMeszaros & ReesI l201lh . Furthermore, the 
conversion is not efficient in uncon fined, time -stationary 
outflows in ideal MHD (e. g., Goldreich & Julian 197d 
iBogovalov & Tsinganos|[l999h . Therefore, models for full 
jet acceleration employing a combination of differential col- 
limation, time variability, or violations of ideal MHD have 
been studied. These result in F, th at increases roughly 
as ex r'^/^ (e.g. iDrenkhahnI HOO?, 'Granot et al.' '2011'; 
iMetzger efaLlbOl lb or ex r^/^ (e.g.. ,McKinnev & Uzdenskv. 



2012h. For illustration, we adopt the calculations of 
McKinnev & Uzdenskv! l l2012h who find the jet dynamics 



well fit by Tj ex r^/^ and E ex r^/^ during jet propagation 
in the star (see their Figure 5). In this case the jet head pres- 
sure falls as ex r^'^/^, so cocoon formation is again inevitable, 
occurring in the range 10^-10^ cm. Note that our adopted 
magnetic -dominated jet has properties — such as particle den- 
sity, energy density, bulk acceleration, and so on — in between 
our conical and funnel jets. 

2.3.2. Survival of Nuclei in theCcocoon 

The freshly formed cocoon consists of an inner region com- 
posed of shocked jet material and an outer region composed of 
shocked stellar material. The two are separated by a contact 
discontinuity. The contact discontinuity is dynamically un- 
stable and_studi£S^rid_that it remains on the order of seconds 
(e.p.. iMizuta & Alovl 12009 ). After that, the outer and inner 
cocoon material mix. We discuss how stellar nuclei may sur- 
vive until this last phase of cocoon evolution. Note that the 
discussions below will also apply to magnetic outflows if they 
interact with the progenitor to produce a cocoon. 

Nuclei in the star ahead of the jet are first shocked by the 
forward shock approaching with velocity ^ /3/i. Spallation 
occurs when the energy of a collision between nuclei and jet 
head protons exceeds the nuclear binding energy « 10 MeV 
in the nuclei rest frame. Therefore, external nuclei survive 
if the jet head velocity is less than /3sp ~ 0-14 in the stellar 
frame. The jet head bulk velocity for a conical (^ = 2) jet. 



«(2) ^nni rl/2 „l/2 
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typically exceeds Psp for radii close to 10^^ cm. For more 
collimated jets, /3sp is reached at smaller radii. For example, 
the jet head velocity of a f = 1 jet of the same jet parameters 
is Ph'^'''' ^ 0.1Lkc,50^^^^9''o,7^^^^- Also note that Ph grows 
faster with radius and the saturation radius is larger There- 
fore, this channel for nuclei loading is limited to jet head radii 
less than 10^-10^^ cm depending on jet parameters. 

Nuclei can also enter the cocoon directly through its bound- 
ary with stellar matter. The cocoon is overpressured and ex- 
pands into the star at velocity /3c given by pressure equilib- 
rium, p^c^Pc^ = Pc, where Pc — Ec/{3Vc) is the radiation 
dominated cocoon pressure and Ec — Lj^ejt — r/c) is the to- 
tal energy dep osited in the cocoon JBegelman & Cioffill 19891 
lMatzneiJl2003b . We approximate the cocoon as a cylinder of 
height r and base length Xc ~ cfict, so that its volume is 
Vc — TTc^Pc^t^r. The cocoon expansion velocity is then 
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which is significantly slower than the head velocity and de- 
pends weakly on jet parameters. For example, even for a 
^ = 1 jet with our nominal jet parameters (i.e., L^c = 
10^°crgs^^, 77 = 100, and tq = 10^ cm), the cocoon ex- 
pands at only /3c ^ 0.04 rio^^^. Furthermore, the flow 
rate of material into the cocoon through the cocoon bound- 
ary can be substantial. The flow rate scales as the product of 
the expansion velocity and surface area. For the jet head this 



is ~ TTXc^cPh- For the cocoon boundary this is ^ 2TTXcrc/3c 
which can be rewritten ~ 2'KXc^cl3h using the approximations 
r K, cjSht and Xc ~ cPct- Thus, the flow rates through the jet 
head and cocoon boundary are Hkely at least comparable. 

Once in the outer cocoon, nuclei must survive photodisin- 
tegration. The photon temperature in the bulk acceleration 
phase is 
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and for a conical (^ = 2) jet this yields 
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where for generality we adopt a parameter Ce ^ 1 to charac- 
terize the fraction of the jet head internal energy that is radi- 
ated by electrons. Throughout, we adopt Ce — 0.1 and show 
dependencies on eg. We quote results for the bulk accelera- 
tion phase because the GRB jet is typically in the acceleration 
phase when it advances through the stellar core, unless it is 
a baryon-rich jet with a small saturation radius. Note that 
for a f = 1 jet of the same parameters, T^^' ^ 160 keV. 
Thus temperatures are insufficient for rapid photodisintegra- 
tion. Also, the cocoon would be turbulent, with flow speeds 
that may reach the sound velocity ^ c/VS. Nuclei being car- 
ried by such flows will exceed the spallation energy threshold. 
As we detail in Section l2.3.3l only relativistic nuclei above a 
critical energy are spalled, because lower energy nuclei are 
expected to thermalize rapidly. 

Next, we must check that nuclei survive potential colli- 
sions with protons in the inner cocoon. In particular, jet pro- 
tons may acquire relativistic random velocities at the reverse 
shock. The protons could be trapped in the fireball by strong 
magnetic fields that are either generated in situ or advected 
from the central engine, and at the reverse shock their direc- 
tions could be isotropized so that they maintain their velocities 
comparable to Vj in magnitude (e.g., see discussions in llokal 
I2OIOI) . In this case, the relativistic protons can cause spalla- 
tion reactions. 

However, the relativistic protons typically lose energy very 
rapidly. For example, relativistic protons lose energy by tt- 
production at an energy loss rate 
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where the quantity in brackets is the jet head particle den- 
sity and ctq is again 3 x 10^^^ cm^. For our canonical jet. 
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for ^ = 2 and ^ = 1, respectively. 7r-production reduces the 
proton kinetic energy down to around 70 MeV. Protons also 
lose energy by processes such as e^-production and Coulomb 
interactions which are particularly important to reduce the 
proton kinetic energy further Since electrons rapidly lose 
their energy by Compton scattering and Bremsstrahlung emis - 
sion, the electrons thermalize and their temperature is lower 
than protons. Thus, ^ 70 MeV protons lose energy to elec- 
trons at a rate (see also Section |2.3.3t 



where In A ^ 10 is the Coulomb logarithm, Ve is the electron 
velocitjQ and n^ — Up if the only electrons in the jet are those 
associated with the ions; if e^ pair production in the jet head 
can increase the proton cooling rate (see Section l372t . For our 



canonical GRB jet, this yields z/. 
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and 



Vpe^^^ ^ 7 X 10^° rg^i s^i for ^ = 2 and C = 1, respectively. 

The short cooling timescales imply that relativistic protons 
occupy only a very thin region around the reverse shock. The 
proton cooling time scale is also much shorter than the life- 
time of the contact discontinuity separating the outer and in- 
ner cocoons. Thus, relativistic protons would have lost much 
of their energy by the time nuclei in the outer cocoon come 
into contact with jet protons in the inner cocoon. Note that 
this conclusion can break down for jets with extremely large 
7], where the optical depth to pp collisions c an be below unity 
and protons will not thermalize (Ilokall2010l) . 

Finally, nuclei must survive in the cocoon. The cocoon tem- 
perature is determined by the energy density of the cocoon. 
We estimate the cocoon temperature as 
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where e,. = 0.1 as before gives tJ^"^' ^ 80 rg ^/^^ eV and 
Tc'^^ - lOOrg-1/2 keV for ^ = 2 and <^ = 1, respectively. 
These are too low for significant photodisintegration. 

In conclusion, external nuclei enter the cocoon through the 
jet head and the boundary between the cocoon and the star 
Those entering through the jet head survive provided the jet 
head velocity is slower than /3sp. This is usually satisfied for 
r < 10^^ cm, but for more penetrating jets can be as limited 
as r < 10^ cm. Those entering through the cocoon boundary 
survive spallation. Once in the cocoon, nuclei survive both 
spallation and photodisintegration. 

2.3.3. Survival of Nuclei During Entrainment 

A classical HL GRB jet is likely to be highly relativistic 
at its central cross section, moving with bulk Lorentz factor 
Tj. Surrounding this is a transition layer where the velocity 
decreases from relativistic to non-relativistic values. Outside 
the transition layer lies the non-relativistic nuclei-rich cocoon. 
The growth time of shear-driven instabilities at the transition 
layer is much smaller than the duration of the jet, and the 
transition layer li kely contains rapi d fluctuations in thermody- 
namic quantities dAlov et al.ll2002l) . Here, we discuss whether 
nuclei are destroyed when nuclei cross into the jet through 
such transition layers. 

When a nucleus in the cocoon moves into the jet plasma, it 
has an extremely short thermal relaxation time corresponding 
to ^ l/{vAp + VAe), where VAp and vac are the energy loss 
rates on jet protons and electrons, respectively. These can be 
written as j/ = 2vs — iy± — t^||, where i^s^ '^±, and j/|| are 
the slowing down rate (or momentum loss rate), the pitch- 
angle diffusion rate, and the parallel velocity diffusion rate, 
respectively. The pitch angle and parallel velocity diffusion 

^ We assume that Te is equivalent to the radiation temperature T/j be- 
cause of Compton scattering. When the temperature is T^ < 10^ eV, 
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refer to the perpendicular and parallel velocity components 
of test particles spreading in velocity space through multiple 
Coulomb scatterings. For sufficiently energetic test particles, 
i/^ and I'll are smaller than vs', in other words, the overall 
slowing down of a test particle is more significant than its 
diffusion in velocity space. We may thus approximate v w 
2i's- The slowing down rate is defined I's — ~(^^||)/'^' 
where v is the particle v elocity and z; || is in the direction of 
the nucleus motion (e.g..' SpitzeiJll956r) . and 

(1 + jAmA/me)ADG{vA/Ve) ,,_, 

^S,Ae = 5 , (17) 

where va is the nuclei velocity and ja is its Lorentz factor, 
Ajj and G{y) are 
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and $(y) is the usual error function. This is the energy loss 
rate on jet electrons; the same expression with e replaced 
by p applies for energy loss on jet protons. The expres- 
sion becomes inaccurate when u^/we ^ InA^e, because 
terms ignored in its derivation become important (" Spitzerl 
Il956l) . However, in our case the jet electron velocity is al- 
ready close to c and this is not a serious concern. Indeed, 
for large y we can approximate G{y) — > l/{2y'^) and we 
find that v agrees numerically with the well-documented en- 
ergy loss rate of relativistic cosmic rays propagating through 
fully ionized plasmas (see, e.g.. Equation (5.3.40) in Sec- 
tion 5.3.8.1 of lSchlickeise3l20()2h . Note that for va < Ve, 
G{y) « {2y)/{3y/Tr) and one obtains Equation (fTsl l. For ref- 
erence, the expressions for the pitc h-angle diffus ion and the 
parallel velocity diffusion rates are (ISpitzerlll956l) 
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These affect the energy loss rates at small nuclei energies. 

In Figure|3] we show the energy loss rate u — 2vs -~v^~v\\ 
for a test Fe nuclei in the jet plasma at r = 10^ cm. The jet ion 
component is assumed to be proton dominated (Zi = 1) and 
the electron and proton temperatures are taken to be equiva- 
lent to the jet radiation temperature; for the bulk acceleration 
phase this is Tj ~ To(r/ro)~'^''^. For high-energy nuclei, en- 
ergy loss on electrons is more important than energy loss on 
protons because of the faster electron velocities. We compare 
the energy loss rates to the spallation rate Vsp — npagpV-pc- 
We do not show the spallation energy loss rate because we 
wish to remain conservative and assume that a single spalla- 
tion event affects the composition. We see that at relativis- 
tic nuclei energies, spallation dominates over Coulomb cool- 
ing rates. However, below -Epc.crit ^ 20 GeV, Fe nuclei 
lose energy before they are spalled, even though spallation 
is energetically allowed. These nuclei lose energy to elec- 
trons on an exponential time scale E'pe oc e"''-*"""**, where 
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Fig. 3. — Energy loss rates and spallation rate for an energetic test 
Fe nuclei moving through the jet plasma. Jet electrons and protons are 
assumed to have temperatures equivalent to the jet radiation temperature. Jet 
parameters are tq = lO'^ cm, L]^^ = 10^" ergs~^, r] = 100, and r = 10® 
cm. Pure proton jet composition with rie = Up has been assumed. Below 
^Fo,crit ~ 20 GeV, Fe nuclei lose energy before they are spalled, even 
though spallation is energetically allowed. 



VAe.crit — l^Ae{EFc,ciit)- 

Importantly, the critical energy -Epc.crit depends very 
weakly on the radius and on GRB parameters. This is be- 
cause in the limit that va > Ve, the Coulomb cooling rate can 
be approximated as 
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for Zi ~ 1. Since the spallation rate is also v^-p oc n^, 
the GRB dependencies largely cancel, and i^Pccrit is always 
0(10) GeV or higheo Equating Vs-p to uac and using the 
general relation n^ = ZiUi yields the approximation 
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(23) 

which is close to the value in Figure [3] Note that this is a 
general critical energy derived from the competition between 
spallation and thermalization in a plasma consisting of elec- 
trons and ions of charge Zi. It can also be applied to, e.g., the 
cocoon (Section r2.3.2b . where Zi > 1. 

Let us denote the bulk velocity gradient at the cocoon-jet 
boundary by df3/dx > 0, where x is in the transverse di- 
rection with the origin at the boundary. As a nuclei moves 
into the jet in the x-direction, collisions with jet protons be- 
come increasingly more energetic, until spallation becomes 
energetically possible. However, the nuclei will also tend to 
thermalize with the jet plasma. Once thermalized, the colli- 

^ If e^ pair production is important, cooling becomes more rapid and the 
ciitical energy becomes higher; see Section 1X2] 
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Fig. 4. — Upper limits on the dfi/dx that are required for Fe nuclei 
that move into the jet medium to survive spallation. Shown for various 
combinations of Ly^^ and r;. For high dj3/dx, the nuclei-jet proton collisions 
become more energetic and nuclei are destroyed. Shading shows this 
spallation region, drawn for our nominal L^^ = 10^'' erg s~^ and 77 = 100. 

sion between nuclei and jet ions will lack the energy to cause 
spallation. Thus, if spallation is slow enough compared to the 
thermalization rate, nuclei will survive. Since the spallation 
rate depends on the collision energy, an upper limit on df3/dx 
can be derived. 

Consider a test nuclei moving in the x-direction with an 
initial velocity va,c ~ 6 x 10^ Tc^5^^'^ cm s^^ corresponding 
to Tc, and let us work in the frame of the jet immediately 
surrounding the test nuclei. After moving a distance Aa;, if 
the nuclei have not yet thermalized with the surrounding jet 
plasma, their velocity perpendicular to x is Lorentz boosted 
by a frame change of c{df3/dx)Ax. As we showed in Figure 
[3] nuclei survival requires the nuclei kinetic energy to be less 
than E'pc.crit- This yields 
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where we have neglected the initial thermal kinetic energy of 
the nuclei as it is much smaller than ii'Fo.crit- Since £'Fc,ciit 
depends weakly on the GRB parameters. Equation (|24] | does 
not show strong parameter dependencies either Now, the dis- 
tance Ax is governed by thermalization. Since vac falls with 
energy, the thermalization distance grows with energy. Larger 
thermalization distances result in larger frame change boosts 
and place stronger constraints on d(3/dx. So we consider the 
thermalization distance at E'Fc.crit, 
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where we have assumed that the energetic nuclei move in a 
straight line (in the x direction) with velocity va.c- This is 
an upper limit, unless nuclei are carried by e.g., bulk flows 
faster than va,c- From Equations ( |24] | and (l25T l we can derive 



an upper limit on the dfi/dx required for nuclei survival. In 
other words, if the velocity gradient is steeper, the collision 
between nuclei and jet ions will become sufficiently energetic 
in the course of the nuclei moving a thermalization distance 
and spallation will occur 

In Figure |4] we show the upper limits on d(3/dx for vari- 
ous GRB jet parameters all for a ^ = 2 jet. The limits be- 
come more stringent with increasing radius, jet Lorentz fac- 
tor, and inversely as the jet kinetic luminosity. These make 
sense, since these provide smaller target densities for nuclei 
to thermalize with and thus must yield stricter upper limits. 
For the same reason, the upper limits for a ^ == 1 jet are re- 
laxed compared to the £, = 2 jet. The breaks in Figure |4]are 
due to the changing dependency of the jet bulk Lorentz factor 
on radius before and after saturation. 

We have made several simplifying assumptions that are 
conservative in nature. First, we have adopted a fixed par- 
ticle density in the entire the cocoon-jet boundary equal to 
the jet particle density. In reality, the density in the cocoon- 
jet boundary will be larger due to a particle density gradient 
from the more dense non-relativistic cocoon to the more tenu- 
ous jet. By fixing to the jet density we have been conservative 
in estimating the upper limit on dj3/dx. We also neglected e^ 
pairs which at small radii will dominate over the e^ associated 
with jet ions. As we detail in Section fTl7\ the production of 



e* pairs aids in nuclei survival because the pairs provide addi- 
tional targets to which nuclei lose energy to, helping nuclei to 
cool on a shorter time scale. This would allow steeper veloc- 
ity gradients. Finally, we have adopted the criterion that the 
spallation rate is smaller This is perhaps conservative in light 
of the fact that a single spallation can still retain intermediate- 
mass nuclei composition. An alternate comparison we could 
have made would be to the spallation energy loss rate. Au- 
thors have adopted other criteria, e.g., Metzger et al.l (tZOl lb 
consider an optical depth < 10 to be acceptable. These relax- 
ations will allow larger velocity gradients. 

We have also neglected magnetic fields at the jet-cocoon 
boundary. Magnetic fields remain highly uncertain but may 
be advect ed from the central engine or amplified in situ. For 
example, IZhang et all (I2OO 9I) find that shear-driven instabil- 
ities induce macroscopic turbulence that leads to the ampli- 
fication of magnetic fields which saturate at an energy frac- 
tion of es ~ 10"'^. If such processes also act efficiently 
inside the progenitor, the magnetic field can be as strong as 
{SneBUhy/^ -- 10^ G at 10^ cm. Strong magnetic fields 
can affect nuclei propagation in different ways depending on 
the field strength and power spectrum, as well as turbulence 
properties. In one limit, nuclei may propagate diffusively in- 
stead of rectilinearly. The effective transverse velocity of nu- 
clei moving into the jet is then reduced, giving more time for 
thermal relaxation and allowing nuclei to survive larger ve- 
locity gradients. In the other limit, magnetic fields may trap 
nuclei. In this case, one may consider nuclei being carried by 
turbulence flows. Such flows can have a variety of velocities, 
lifetimes, and mixing efficiency. They may reach up to the 
sound speed ^ cj\fi. From Figure |3] nuclei moving at speed 
c/\/3 is in the spallation regime. However, it is marginal, and 
nuclei may still survive considering the conservative assump- 



tions described above. 

Finally, we have implicidy assumed that baryon entrain- 
ment occurs efficiently while the jet propagates through the 
star. We discuss the implications of the uncertain entrainment 
efficiency in the discussions section. 

2.3.4. Expected Velocity Gradient at the Jet-Cocoon Boundary 

The dynamics of relativistic GRB jets propagating through 
their progenitor stars have been extensively studied both 
analy t ically and numerically (Aloy et al . 2000; Zhang et aU 
2003t iMatzneJ 120031: IZhang et al.1 120041: iMizuta et al.ll200a 



Morsonv et al.l2007l) . While simulations confirm that the bulk 



velocity drops sharply at the jet boundary, with Vj decreasing 
to order unity, a meaningful comparison with Figure 2]is not 
yet possible. Even with state-of-the art numerical simulations, 
the desired transverse resolutions have not yet been reached. 
Fo r example, in the s imulations of relativistic jet propagation 
bv lZhang et al.l (l2003l) . the angular resolution is 0.25 degrees, 
i.e., a few 10^ cm at radii of ~ 10^ cm. Finer features are 
expected an d seen with increasing simulation resolution (e.g., 
iMizuta etal. 2006) . but simulations do not yet reach compa- 
rable scales as Figured 

However, it is clear that the velocity gradient should be con- 
nected to instabilities at the boundary. Alov et al. (2002) stud- 
ied the stability properties of GRB jets propagating through a 
dense medium like the progenitor star. Through both numer- 
ical simulations and linear stability analysis, the authors find 
that the shear at the jet-cocoon boundary is responsible for in- 
stabilities that grow on very rapid time-scales. They show that 
the growth rate of non-homogeneous radial perturbations in- 
creases with the velocity gradient, i.e., vku ^ Vj\V'\, where 
\V'\ is the initial velocity gradient. Therefore, a large \V'\ 
would cause the rapid growth of instabilities which would 
tend to smooth the velocity gradient to a smaller final veloc- 
ity gradient df3/dx. Recall that nuclei survival requires the 
velocity gradient d/3/dx to be small over distance scales of 
^ Axthcrm- Thus, having a large initial gradient \V'\ may 
even help provide the necessary shallow d/3/dx for nuclei sur- 
vival. For example, adopting the upper limit ^ 0.2 cm^^ 
(at 10® cm; Figure |4]i for \V'\, the instability growth rate 
z^KH ^ 6 X 10^^ s~^ is much faster than the spallation and 
energy loss rates, so it is likely instabilities would grow and 
reduce the velocity gradient. For much smaller initial velocity 
gradients, instabilities may not have sufficient time to grow, 
but the small initial velocity gradient would anyways work 
positively for nuclei survival. 

3. EFFECTS OF BULK ACCELERATION AND DISSIPATION 

In the previous section, several candidates of the origin of 
jet nuclei were identified. However, even if nuclei are suc- 
cessfully loaded into the jet, they may be destroyed during 
the evolution of the jet because of bulk acceleration and/or 
dissipation. Here we discuss some of these effects. 

3. L Neutrons during Bulk Acceleration 

Initially, the neutron component of the jet is well coupled to 
the ion component by elastic collisions with a small relative 
velocity /3 ^ 1 that is insufficient to cause nuclei to break up. 



However, the neutrons will lag behind the ions in the plasma 
during the bulk acceleration phase, if the collisions cannot 
keep up with the jet expansion. The velocity lag of neutrons 
relative to charged ions is AF/F ^ Tcoii/roxp, where r„ = 
r — AF is the neutron Lorentz factor, TcoIi — l/("-iO'i/3c) 
is the comoving collision time-scale, ai ~ <7oA'^/^/(3 with 

r/{Tjc) is the comoving 
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expansion time-scale. The neutron-ion relative velocity is 



p 



(3-Pn 
1 - I3(3n 



r„ 



r 



r 



1.^70^2/3 



(26) 



where /3„ is the neutron velocity. The relative velocity in- 
creases with radius during the bulk acceleration phase as 
(3 ex L]ic^^r^^~^7]. Therefore, we require that at the end of 
the bulk acceleration phase (r ~ rsat), the relative velocity is 
smaller than figp- This yields the condition for no spallation 
of 
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for ^ = 2. To be conservative we have assumed an Fe com- 
position for the jet: this leads to a larger /3 than a proton jet, 
because of the longer collisional timescale. For a ^ = 1 jet, 
the saturation radius is larger by a factor 7/ and the expansion 
time-scale is longer Also, the higher jet density results in a 
faster collision time-scale. Both serve to relax the nuclei sur- 
vival condition, yielding 77 < 13G0iko,50^ '^o,?^^^'^- Similar 
conclusions hold for our adopted magnetically dominated jet 
model: the survival condition is 77 < 480Lke,5o^ ' o,7~^^'^- 
If the above conditions are not satisfied, the neutron density 
is typically high en ough that nuclei are spalled very rapidly 
(lBeloborodovll2003l) . We plot the constraints for a ^ = 2 jet 
in Figure |5] 

An alternate way for nuclei survival is if neutrons decou- 
ple from the accelerating plasma. We define the decoupling 
radius as the jet radius when the separation between neutrons 
and ions grows larger than the radial width of the ejecta, cAT, 
where AT is the jet pulse duration. In the lab frame, the sep- 
aration is 
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(for ^ 7^ n/5) which implies that decoupling is more likely to 
occur in low Lkc and high 77 jets; for example, it does not oc- 
cur for our nominal ike = 10^°ergs^^ and 77 = 100 jet. Let 
us therefore consider a jet of luminosity Lkc = lO^'^ergs"^ 
and terminal Lorentz factor rj = 100. Note that from Equa- 
tion (l27T i. it can been seen that spallation will occur in such a 
jet. Now, the spallation radius (defined as when (3 = (3sp) and 
the decoupling radius are, respectively. 
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both for C = 2. Requiring that decoupling precedes spallation 
yields 
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Fig. 5. — Constraints on jet parameters for Fe nuclei survival during jet 
bulk acceleration, including spallation caused by collisions with neutrons 
(lines increasing with r\) and effects of neutron decoupling (lines decreasing 
with 1], labeled), for several values of tq. Shaded region results in nuclei 
spallation (shade shown for our nominal ro = 10^ cm and AT = 1 s). 



As expected, decoupling becomes important for sub-luminous 
jets. For example, a small Lkc strongly limits the range of r\ 
for which spallation does not occur [Eq. (|27tl, but the con- 
sideration of neutron decoupling allows a much wider range 
to be acceptable [Equation dSTTil. However, we caution the 
strong dependencies of Equation ( [3T] i on GRB parameters. 
We conclude that decoupling may play a role but its impor- 
tance is strongly parameter dependent. Note that for a ^ = 1 
jet the particle density is larger and the relative velocity /3 is 
significantly slower, so that decoupling typically does not oc- 
cur while the jet is propagating through the stellar core. In 
Figure|5] we show constraints on ike and r\ for several values 
of ro, for a ^ = 2 jet with AT = 1 s. 

In conclusion, nuclei can survive over a significant range of 
Lkc and r\ because of slow neutron-ion relative velocities. For 
low Lkc (and AT), decoupling may occur before spallation so 
that a wider range of r] is allowed. For canonical values of ike 
and 77, collisions just reach spallation thresholds. For more 
collimatedjets, a larger GRB parameter range is allowed. 

3.2. Dissipation and Particle Heating at Oblique Shocks 

Until now we have parameterized the jet cross section by 
Equation (JUi which can be adapted to any morphology, e.g., 
conical (constant jet opening angle) or funnel (jet opening an- 
gle ex r^^'^) behavior. However, recent numerical and an- 
alytic works of jet propagation show that the jet morphol- 
ogy is not a free parameter but correlated with the jet's in- 
te ractions with the cocoon. According to the analytic work 
of iBromberg et al. ( 201 lb. which is ba sed o n the jet-cocoo n 
solutions of BegelmajL&Cioffil (Il989l) and iMatzned (l2003h . 
the jet morphology is determined by the jet energy density, 
its opening angle, and the density of the external medium. 
For a jet propagating through a dense material such as the 
progenitor star, the cocoon pressure is sufficiently strong and 



compresses the jet. Oblique collimation shocks form which 
generate the pressure needed to counterbalance the cocoon's 
pressure. The shocks can form early in jet propagation but 
converge as the entire jet is shocked at a radius of approx- 
imately rsh ^ 0-lfft., where r^, is the jet head radius (e.g.. 
Figure 2 of'Bromber g et al.ll201 ll) . 

For a sufficiently fast jet (Tj > l/9j), th e Lorentz factor 
of the jet after being shocked is F^ ^ l/9j (IBromberg et al.l 
l201lh . which for our GRB parameters Oj « 0.13f2o,-i^^^ 
and thus Fj ~ 8ilo,-i • The jet ions passing the shock 
may obtain high random Lorentz factors corresponding to the 
relative Lorentz factor ^ rj/{2rs) in a similar way to the 
reverse shock (Section [2. 3. 2t . Such ions would collide with 
other jet ions and any nuclei could be spalled. Here, we con- 
sider the spallation rate and energy loss rate to see whether 
jet nuclei are actually destroyed. To place conservative limits, 
we consider the fate of Fe ions within a jet plasma consisting 
of mainly protons. This assumption yields the highest spalla- 
tion target density; the spallation rate is then i^sp — npa^pC, 
yielding z^sp^^) ~ 3 x 10^ikc,5or-sh,9"^?/2"^ s^^ (( = 2) or 



-3/2, 



To estimate the Fe-electron collisional cooling rate, we con- 
sider in addition to the electrons associated with jet protons 
the production of e^ pairs. This directly increases the Fe cool- 
ing rate since i^Ae oc n^ . In fact, it is important for the survival 
of relativistic nuclei, since we showed in Section 12.3.31 that 
only mildly relativistic Fe ions can survive when electrons 
associated with jet protons are considered. For a radiation 
tempe rature Tr the equilibrium pair density is (IShemi & PiranI 
1 19901) 



n± 



4.4xl030(r,/me)3/2g-me/T._ 



(32) 



We estimate the pair density in the jet by substituting the jet 
radiation temperature. For example, for a ^ = 2 jet in the bulk 



acceleration phase, T^ 



(2) 



180ike,50^^*r9"^ keV and the 



resulting pair density exceeds those associated with jet pro- 
tons by a factor of approximately 10^ at 10^ cm. However, by 
^ 10^^ cm the Tj has fallen too low and pairs make a negli- 
gible contribution to the total electron density. For a narrower 
jet, the jet temperature remains high longer so that pairs make 
meaningful contributions out to larger radii, e.g., for ^ = 1 
and canonical jet parameters, out to ~ 10^*^ cm. Note that the 
jet radiation temperature can be high but photodisintegration 
is still slower than the spallation rate. 

Let us first discuss when the fast jet approximation of 
iBromberg et al.l (|201 lb is valid, i.e., Tgh is larger than 8 x 10^ 
cm (^ = 2) or 6.4 x 10^ cm (^ = 1), both for our canonical 
choices of ro and ^Iq. Whether collisional cooling prevents 
spallation depends quantitatively on the Lorentz factor of the 
Fe ions and the electron density. For example, at Vgh — 10* 
cm, the Fe ion cooling rate is I'Ae'"^^ ^ 1 x 10^^£'fc.12^^ s~^ 
for our canonical GRB parameters and £, = 2. Here we take 
the electron temperature to be equivalent to the radiation tem- 
perature, because electrons reach local thermodynamic equi- 
librium on much shorter time scales compared to the ion- 
electron energy loss time scale. Comparing i^^ie*^^'' to the spal- 
lation rate i^sp'^^^ we see that Fe ions will undergo spallation 
if their energies are above Epe ~ 4 x 10^"^ eV, or a Lorentz 
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Fig. 6. — Constraints on jet luminosity forFe nuclei survival in collimation 
shocks. Shaded region results in nuclei spallation (shown for a ^ = 2 and 
r] = 100 jet). A narrower jet is more favorable for nuclei survival due to 
its higher temperature (dashed line). The r; = 5 jet (dotted) is simply an 
illustration where the post-shocked jet Lorentz factor is deliberately fixed so 
that nuclei are spalled; it is possible that the post-shocked jet is faster in which 
case nuclei will not be spalled. For all other jets, the post- shock dynamics 
are ca lculated according to the collimation shock model of iBromberg et al.1 
i2011h . see the text. 

factor of ~ 800. Thus for our canonical parameters, Fe ions 
easily survive at Tgh = 10® cm. 

However, as rsh increases, the pair density drops and sur- 
vival becomes increasingly difficult. In Figure |6] we show 
the jet luminosity required for nuclei survival as a function of 
Tsh- This is determined by calculating the Fe Lorentz factor as 
a function of rsh, and then the necessary pair density so that 
the energy loss rate i^Ae is faster than the spallation rate I'sp- 
It is clear that higher luminosities are required for survival at 
larger radii, because the temperatures are lower For a similar 
reason, a narrower jet is more favorable. A larger value of rj is 
only slightly more favorable mainly due to the smaller spalla- 
tion target density it implies. Note that the importance of pairs 
can be noticed: for Liad = 10^°ergs^^, nuclei survive out to 
Tsh « 10®-^ cm (C = 2) or r^h ~ 10^" cm (^ = 1); these radii 
are the radii where pairs cease to make a contribution. At the 
small radius end, the required luminosity suddenly drops be- 
cause the Fe Lorentz factor tends to unity as the radius where 
the fast jet appro ximation become s invali d is reached. 

The model of IBromberg et al.l (1201 ll) cannot be applied 
for radii where the fast jet approximation is not valid, i.e., 
r^h < 8 X 10^ cm (^ = 2) or r^h < 6.4 x 10* cm (^ = 1), 
for our parameters. However, that is not to say oblique shocks 
will not occur in this region. The stellar pressure is typically 
much greater than the jet pressure so that the jet is pressure 
confined by the star, a setup conducive to oblique shock for- 
mation. In fact, much of the early jet phase can be shocked 
and supply the cocoon material. We do not attempt to model 
the dynamics of these jets. However, we can reason that due 
to the high jet radiation temperatures, it is unlikely nuclei are 
destroyed. In these radius ranges, the pair density exceeds the 



proton density by ^ 10'' and cooling is rapid. For example, 
at r^h = lO^'^ cm (^ = 2) or rsh = 10® cm (^ = 1), the 
cooling rates are vac ~ 5 x 10^'^i?Fc,i2^^ s^^, and spallation 
only occur for highly relativistic nuclei with Lorentz factors 
greater than ^ 1000. 

Similarly, jets with rj < l/6j do not satisfy the fast jet ap- 
proximation, but oblique shocks are nonetheless expected as 
the jet expands into the stellar material and later on into the 
cocoon material. In this case we treat Tg as a free parame- 
ter If Ts ^ 7], then nuclei are not spalled because the relative 
Lorentz factor between the jet and post-shocked jet is small. 
On the other hand, if Tg ^ rj, the nuclei can attain mildly 
relativistic energies and can be spalled. For example, nuclei 
will survive in a 77 = 5 jet if Tg ~ 2.3 or higher In Figure 
|6]we shown as an example the case where Tg is fixed such 
that nuclei are spalled (shown for Tg = 2; dotted). As ex- 
pected, the smaller relative Lorentz factor results in a weaker 
requirement than our canonical (r/ = 100) jet of the same 
morphology (dashed). 

Oblique shocks can also occur after jet collimation as the 
jet propagates through the progenitor, but the survival of nu- 
clei depends largely on the post-shocked jet Lorentz factor. If 
Tg ^ rj as is usually expected, nuclei are likely not destroyed. 

Note that in all the above estimates we have conservatively 
assumed a proton-dominated jet and considered the survival 
of a (minor) Fe component. If the jet is mostly Fe dominated, 
the spallation target density decreases by ^4 = 56 and im- 
proves Fe survival prospects. Also, we have conservatively 
considered the spallation rate. If we relax this and consider 
instead the spallation energy loss rate, i.e., we allow more 
than one spallation event, the limits for nuclei survival will 
be relaxed. 

To conclude, jet nuclei may be spalled if jet ions obtain 
relativistic random velocities at oblique collimation shocks. 
However, when the shock radius Vgh is small, the pair density 
in the jet can be sufficiently high that nuclei energy loss is 
faster than spallation. For our canonical jet parameters, this 
condition is realized for shock radii less than ^10^° cm for an 
initially narrow jet (which may be realized if the jet is initially 
confined by the rotational funnel geometry of the progenitor) 
or less than ^ 10®'^ cm for an initially conical jet. Nuclei that 
are entrained at larger radii (i.e., r > r^h) may not affected by 
the oblique collimation shocks and are not spalled. Jets with 
small rj imply small Fe ion Lorentz factors and are thus more 
favorable for survival. 

3.3. Dissipation and Particle Acceleration at Emission 
Regions 

GRB emissions typically consist of prompt gamma-ray 
emission and afterglow emission. The former is attributed to 
non-thermal or quasi-thermal emissions produced via some 
internal dissipation, e.g., internal shocks or magnetic recon- 
nections. In the standard optically thin synchrotron scenario 
of the internal shock model, inhomogeneities in the jet cause 
internal shocks within the jet, where the relative bulk kinetic 
energy is dissipated and accelerated electrons radiate gamma 
rays. In this scenario, it is natural to expect that ions are ac- 
celerated as well. The maximum energy to which ions are 
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accelerated is roughly determined by requiring that the accel- 
eration timescale is shorter than the dynamical and any energy 
loss timescales. Previous works have shown that the accel- 
eration of nuclei to UHECR energies is possible for certain 
ranges of GRB parameters w hen internal dissipation happens 
in the optically thin regime (iMurase et al.ll2008t I Wang et aU 
I2OO8). 

The conditions required for the accelerated nuclei to sur- 
vive in the prompt GRB photon field are much more stringent 
than those for acceleration. The comoving photodisintegra- 
tion time scale for a nucle us moving through an isotropic pho- 
ton background is (e.g.. ,Mura"se et al.bOOslilWang et al.l2008l) 



M7 — 



1 c p , , , f^ ^ 1 dn 

7 = ^rr/ deecjA^{e) ^e -^ — , (33) 

^7a "'eth Je/2'YA ^ "^ 

where quantities such as e are defined in the nucleus rest 
frame, 7^ is the nucleus Lorentz factor, Cth w 7.6 MeV is the 
threshold photon energy for photodisintegration of an iron nu- 
cleus, and dn/de is the differential photon spectrum. Adopt- 
ing a comoving broken power law for the prompt GRB photon 
spectrum. 
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(34) 

where £& = £b,obs/'7 is the break energy, the factor 5 in the 
denominator arises because the luminosity at the break energy 
is approximately 1/5 the total luminosity, Eb.obs ~ 1 MeV, 
Emin = 1 eV, and Emax = 10 McV, the comoving photodisin- 
tegration timescale is estimated to be 
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300eg iiL^^Vi5?72eb,obs,6 s, 



(35) 



where we have adopted ^ = 2 since the jet becomes conical 
after jet break out. We require that Tcxp < ^^7 for UHECR 
nuclei to survive photodisintegration in the system, which 
yields the constraint 



ri > 110e^;_ii 



1/2 ^1/2 -1/2 -1/2 



'b,obs,6- 



(36) 



We plot this in Figure Q as a function of the dissipation 
radius, i.e., emission radius. We emphasize that the optical 
depth for photodisintegration highly depends on the dissipa- 
tion radius that is estimated to be ris w 2c St rj^ for the in- 
ternal shock ca se. Typica l internal shock radii are ^ 10^'^- 
lO^^-^ cm (e.g.,|Naki&^Piran 2002), so survival of UHECR 
nuclei is possible only for relatively large dissipati on radii, al- 
thoug h this is relaxed for sub-luminous GRBs ( Murase et al.l 
2008h . Typical radii for magnetic dissipation scenarios simi- 
larly span a wide range but can reach up to ^ 10^^ cm ( e.g.. 



lDrenkhahnl2002HGranot etai]|201 itlMetzger et al.ll20TTl) . In 
Figure I2] we show these ranges of dissipation radii for refer- 
ence, noting that there are substantial uncertainties. Note that 
in the photospheric emission scenario, the emission radius is 
typically expected to be much smaller, where high-energy nu- 
clei break up via both spallation and photodisintegration and 
instead high-energy neutrinos are generated (Mu rase 2008). 

Dissipation and emission also occur later on, when the 
jet expands into the interstellar medium and sweeps up suf- 
ficient baryons to enter the Blandford-McKee phase. The 
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Fig . 7 . — Constraints on jet r; and internal dissipation radius for Fe survival 
at the dissipation region. GRB jets of smaller luminosity and larger tj 
allow smaller internal dissipation radii to satisfy nuclei survival conditions. 
For illustration, we indicate typical dissipation radii for the internal shock 
scenario (rig) and the magnetic reconnection scenario (rmag) by arrows. 

forward and reverse shocks are collectively referred to as 
external shocks, where ions and electrons can be acceler- 
ated. The observed non-thermal photons are attributed to 
synchrotron emission from relativistic electrons accelerated 
at such shocks. Compared to the case of the prompt emission, 
nuclei interact with softer photons, and it has bee n shown that 
the survival of UHECR nuclei is also possible JMurase et al.l 
l2Q08l:IWangeraDl2008h . 

4. DISCUSSIONS AND CONCLUSIONS 

In this work, we focus on the fate of nuclei in a relativis- 
tic jet accompanying a CCSN. We consider three sources of 
jet nuclei, (1) loading at the jet base during jet launch, (2) in 
situ explosive nucleosynthesis, and (3) entrainment of exter- 
nal nuclei during jet propagation. We first discuss the condi- 
tions for nuclei survival in each of them. For initial loading, 
we find that nuclei survive if the jet launch radius is greater 
than ro ^ 10® cm or the initial radiation luminosity is less 
than Ljad.o ^ 10^® ergs~^ (Figure|2]l. These conditions may 
be satisfied by magnetic models of classic GRBs or models 
of sub-luminous GRBs. If nuclei are destroyed at jet launch 
into free nucleons, the material is fused into nuclei, but the 
freezeout abundance contains significant heavy nuclei only 
for low entropy jet s (IBeloborodovl l2003t llnoue et aP 120031: 
iMetzger et al.ll201 lb . Finally, whether external nuclei survive 
during entrainment into the jet depends critically on the ve- 
locity gradient at the cocoon-jet boundary (Figure 2)- If the 
gradient is too steep, collisions between external nuclei and 
jet ions become too energetic and external nuclei are spalled. 
On the other hand, if the gradient is shallow enough, the nu- 
clei thermalize with the jet plasma before spallation and sur- 
vive. Since the growth rate of shear-driven instabilities at the 
cocoon-jet boundary is larger than the nuclei energy loss or 
spallation times, the velocity gradient is expected to be shal- 
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low enough to allow nuclei to survive. 

Next we investigate the conditions for nuclei survival once 
loaded into the jet. First, collisions with neutrons become en- 
ergetic during the jet bulk acceleration phase and can cause 
nuclei spallation if the GRB terminal Lorentz factor is larger 
than r/ ^ 220Lke,5o^ ^"0,7"^^^ (for a conical jet) or 7] ^ 
1300Lkc.5o^ ''o,?"^^^ (for a funnel jet); for smaller 77, spalla- 
tion is energetically prohibited. Also, we identify the param- 
eter space where the neutrons decouple before causing spalla- 
tion. We find that it is strongly parameter dependent but tends 
to help with small Lke jets (Figure |5]l. Second, we show that 
nuclei can survive at oblique shocks if e^ pairs are produced 
such that the nuclei collisional energy loss rate is competi- 
tively high. Such conditions are typically realized for shock 
radii r^h less than ^ lO'^-lO^" cm (Figure |6]l. A smaller 77 
works positively for nuclei survival, as does a post-shocked jet 
Lorentz factor Tg that is close to the pre-shocked jet Lorentz 
factor 

Based on the above results, we can consider the sources and 
survival of nuclei in three distinct sets of GRB models. These 
are summarized below and in Table [T] 

First, we conclude that nuclei from initial loading and ex- 
plosive nucleosynthesis are disfavore d in th e HL GRB fire- 
ball scenario (see also Beloborodovl 120031) : instead, nuclei 
may come from entrainment during jet propagation. Once 
loaded, nuclei may survive collisions with neutrons. Survival 
in oblique coUimation shocks is possible while r^h is less than 
^ 10^-10^" cm. After r.p.m.sh reaches these values, nu- 
clei entrained at smaller radii (r < r^h) may be spalled as 
they cross the shocks, while nuclei that are entrained at larger 
radii (r > rgh) may survive. Successfully entrained nuclei 
can survive the emission region provide d the dissipation ra- 
dius is large eno ugh (Figure [Tj see also iMurase et al.luOOSt 
IWang et al.ll2008h . 

On the other hand, multiple sources of nuclei are possible 
for magnetic-dominated scenarios of HL GRBs. For example, 
nuclei may survive at initial loading, because for a given GRB 
luminosity the initial radiation energy can be lower Even if 
nuclei are spalled at loading, the subsequent explosive nucle- 
osynthesis can lead to Fe-group elements and possibly beyond 
(JMetzger et al. 2011). Entrainment can work too, if as in re- 
cent models the jet is collimated by a cocoon. The slower 
bulk acceleration compared to fireballs means collisions with 
neutrons are generally less destructive, and oblique shock are 
expected to be weaker, although definite statements are de- 
pendent on the specific model of magnetic dissipation. 

Similarly, there are multiple potential nuclei sources for 
sub-luminous GRBs and jet-driven CCSN. Their smaller lu- 
minosities and bulk Lorentz factors work positively for nu- 
clei su rvival at initial lo ading and for explosive nucleosyn- 
thesis (llnoue et al.ll2003b . Although the low ike works nega- 
tively for spallation by neutrons, this is largely compensated 
by the lower ry; as a result, nuclei are typically not destroyed 
(for the smallest Lko, neutron decoupling can also be effec- 
tive). The low T] also works positively for nuclei survival in 
oblique collimation shocks, because the nuclei Lorentz fac- 
tors are small. However, note that nuclei can still be spalled 
at large Vgh if ^s -^ V- Finally, though uncertainty is large. 



TABLE 1 
Summary of Nuclei Survivability 





HLGRB 

Fireball 


HLGRB 

Magnetic 


LLGRB 


Source: initial loading 
Source: nucleosynthesis 
Source: entrainment 


N 

N 

maybe 


Y 
Y 

maybe 


Y 

maybe 
maybe 


n collisions 
Collimation shocks" 
Emission region'' 


Y 

maybe 
Y 


Y 

maybe 

Y 


Y 
maybe 
Y 



Note. — "N" denotes not possible for typical parameters; "maybe" pos- 
sible depending on model, and "Y" possible for typical parameters. Typical 
parameters are taken to be Lko = lO^'^ergs"^, rj = 300, ^ = 2, and 
ro = 10^ cm for HL-GRBs and Lkc = lO'*'' ergs"!, r; = 10, g = 2, 
and ro = 10^ cm for LL GRBs. For entrainment, we consider the currently 
unconstrained jet-cocoon velocity gradient to be model-dependent. 
" For all cases, other oblique shocks may lead to nuclei destruction depending 
on shock velocity. 
'' For emission radius 10^^ cm. 

survival at dissipation is easier than in HL G RBs (Figure |7j 
see also lMurase eTaPbOOsHWang et al.ll2008h . 

As seen in this work, jets accompanying CCSN may con- 
tain heavy nuclei that originate from the stellar core, disk 
wind, and/or low entropy (kinetically or magnetically dom- 
inated) jets. A high abundance of nuclei is attractive in view 
of recent reports of a heavy-ion dominated composition of the 
highest-energy UHECR. If AGNs are the sources of UHE- 
CRs, rigidity-dependent acceleration models naturally pre- 
dict a heavy-ion composition at the highest energies, but this 
seems inconsistent with the null observation of an excess of 
protons at lower energies (Abreu et al. 2011b). On the other 
hand, nuclei-rich UHECRs could be realized in GRBs and 
CCSN, as long as the accelerated particles are injected with- 
out being broken up. 

However, GRBs and relativistic CCSN scenarios may also 
have several issues such as an energy crisis and cosmic-ray 
escape problems. That an energy-crisis problem may be se- 
rious for classical HL G RBs has been claimed in view of 
Fermi observations (e.g., lEichler et al.l 12010) . although the 
universality of this claim depends on the uncertain local rate 
of classical GRBs, the efficiency of cosmic ray accelera- 
tion, and the slope of the injected UHECR spectrum ( e.g., 
iLe & DermeJl2007l iMurase et al.ll2008[ IWaxmanll2010l) . It 
is interesting to note that magnetic GRB models may pre- 
dict flatter accelerated cosmic ray spectra which could help 
to avoid an energy crisis (e.g.. iRomanova & Lovelacdl 19921; 
iMetzger et al.ll201 lb . although it remains to be seen whether 
it can also offset the issue of smaller baryonic loading in mag- 
netic jets. Furthermore, sub-lu minous GRBs may contribute 
to the observed UHECR flux dMurase et al.|[2006 ). The ex- 
pected rate of CCSN with relativistic ejecta (whether it is 
jet-like or not) seems co mparable to that of sub-luminous 
GRBs dChakraborti et al.ir201 1 ) . and both jet and relativistic 
CCSN ( which may be jet-driven) can provide UHECR nu- 
clei (e.g..'Muras e et al.ll2008l TWang et al.ll200 8'). 

We emphasize that we have attempted to place conserva- 
tive constraints on nuclei survival. This means we have not 
always used the same jet composition in each investigation. 
For example, when we consider neutron collisions we adopt 
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a nuclei-dominated jet so that the relative velocity is high- 
est and spallation occurs most readily; and when we consider 
oblique shocks we adopt a proton-dominated jet so that the 
spallation rate is maximal. In the same spirit, we assume that 
a single spallation effect must be avoided, when in fact a few 
spallations may still maintain a heavy or intermediate nuclei 
composition depending on the initial composition. Relaxing 
our criterion and allowing, e.g., the spallation optical depth 
to be a few, relaxes the nuclei survival constraints described 
above. This would naturally lead to a final composition of 
more intermediate nuclei. 

We have been simplistic on several issues. We have only 
considered two kinds of jets, and also not considered pro- 
genitor dependencies. But the most important is our implicit 
assumption that loading and entrainment of external baryons 
occur efficiently. Jets accompanying CCSN can have a wide 
range of baryonic content, with GRB jets among the clean- 
est with fewest baryons. GRB jet models achieve this in dif- 
ferent ways. For example, magnetic-dominated jets that are 
powered by energy release on field lines where baryons are 
confined can plausibly conta in very few initial baryons (e.g., 
iLvutikov & Blandfordll2003l) . Thus, although nuclei are more 
likely to be present at the launch sites of magnetic-dominated 
jets (Table [Til, they may not be efficiently loaded. However, 
we stress that such a statement remains highly model depen- 
dent. Later entrainment is attractive in the sense that jets 
must load baryons at some point, and propagation through the 
baryon-rich progenitor is an unavoidable course of jet evolu- 
tion. The entrainment efficiency is expected to depend on the 
instability mechanism, thermodynamic parameters of the co- 
coon and jet, magnetic field, and so on, and thus evolve. To 
make quantitative estimates of these and other effects requires 
a detailed model of entrainment and is beyond the scope of 
the present paper In this work, we investigated the process 
by which nuclei may survive during loading. Although the 
specific GRB model may change in the future, the physical 
processes by which nuclei survive should still hold. 

We can however speculate about the composition assum- 
ing entrainment is the main process of jet baryon loading. In 
an oversimplified picture where entrainment occurs equally 
efficiently during the entire propagation of the jet, the jet 
composition will be similar to the total matter swept up in 
propagation. However, this is likely strongly modified. First, 
explosive nucleosynthesis is expected to increase the nuclei 
abundance of the swept up material. Indeed, high ^^Ni 
masses have been observationally supported by nearby highly 
luminous CCSN dlwamot o et alJll998l: IWooslev et alJ ll999l: 
IWooslev & Hegej2003h . and observational hints of jet-mixed 
nuclei have been observed (e.g.,|Smith et al. 2012). Second, 
nuclei tend to be more easily destroyed during the later phases 
of jet evolution, in particular as the jet head velocity exceeds 
/3sp ^ 0.14, and/or when the oblique shock radius exceeds 
Tsh ^ 10^*^ cm. Nuclei that lie at radii > 10^° cm, mainly 
intermediate nuclei such as carbon, are thus more likely to be 
spalled before entering the cocoon. 

Finally, it is relevant to identify observational signatures 
of heavy nuclei. For protons, neutrinos are produced when 
accelerated protons interact with photons and/or target pro- 



tons and produce charged pions. The neutrino signatures 
have been well studied in various contexts. If collision- 
less shocks can be formed inside the progenitor star, GeV- 
TeV neutrinos might be expected during jet propagation 



(e.g..lMeszaros & Waxmanl200ll:lRazzaaue et alj2003ll2004 



'Ando & BeacomI 120051: iHoriuchi & Andol l2008t llocco et alj 
2008). One may expect TeV-PeV neutrino emission as- 
sociated with photospheric and/or shock breakout emis- 
sions, where the formati on of colhsionless shocks is ex - 
pected (Murase 2008; M urase et all 1201 ll iKatz et alj|201 ih . 
At larger dissipation radii above t he photosphere, PeV-EeV 
neutrinos can be produc e d (e.g., Waxman & Bahcap 1 19971 
iDermer & AtovanI 120031; iMurase & Nagatakil bOOell Al- 
though nuclei also produce neutrinos, the condition that nuclei 
survive limits the amount of neutrino-generating collisions, so 
the neutrino flux is typically small and not so easy to detect 
by IceCube (e.g.. lMurase & Bea com 2010b). We emphasize 
that in heavy nuclei scenarios, the GRB-UHECR hypothesis 
is consistent with the latest upper limits placed by IceCube 
dlceCube Collaboration et aLll2012l) . Also, it would be hard 
to distinguish the n eutrino signature from those of protons 
dMurase et al.ll2008h . 

The existence of nuclei can also be probed by iden- 
tifying atomic or nucle a r line emissions. As in dis- 
cussed iMeszaros & ReesI (120011) . the cocoon material pro- 
duced by the jet would break out of the stellar enve- 
lope, where iron-enriched clumps may be shed by a UV/X- 
ray continuum and an Fe line luminosity of ~ 4 x 
10**^ ergs"^ (Afb/lG"^MQ)a;Fo is expected. Here Mb is the 
bubble mass and xpe is the Fe mass fraction. Alternatively, the 
Fe line emission may be caused by X-rays produced by a con- 
tinuous but decaying jet (Rees & Meszaros 2000) . Though 
the details are uncertain, if the high abundance of heavy nu- 
clei is realized. X-ray line features might be constrained or 
even detected by current and future X-ray satellites such as 
Astro-H. Furthermore, it is possible that boosted atomic pho- 
tons could be expected in the GeV range for accelerated nuclei 
dKusenko & Voloshinll20Tlh . 

Nuclei can also emit ^ MeV gamma rays (in their rest 
frame) through their excitation states. Such nuclear gamma 
rays are not so easy to detect for distant GRBs, but they may 
be detected for nearby and/or luminous events. In particular, 
nuclear gamma rays may be boosted to the TeV energy range 
when nuclei are accelerated. For example, the radioactive iso- 
topes ^^Ni and ^^Co lead to the production of MeV nuclear 
gamma-ray lines after they are excited via, e.g., e lectron cap- 
ture or positron emission (for de tails, see, e.g., iMilne et akl 
I2OO4I: iHoriuchi & BeacomI l20l d). These nuclei can be en- 
trained andJ^adJoGeV-TeVgamma rays in GRBs and hyper- 
novae dioka & MeszarosluOlOl) . The time scale of emission is 
also extended by the nuclei Lorentz factor to ^ lO"" yr 7^.5 
for ^^Co (where 7^ is the Lorentz factor of the nuclei in the 
observer frame), and detection is possible only for Galac- 
tic events. On the other hand, nuclei interact with low- 
energy photons in the source via the photodisintegration pro- 
cess, leading to subsequent (almost prompt) ^ 0.2 TeV ^a,5 
gamma rays from daughter excited nuclei. This signal is use- 
ful as a unique probe of nuclei acceleration as well as syn- 
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chrotron and inverse-Compton ga mma rays from pairs gener- 
ated via the Bethe-Heitler process (iMurase & Beacoml2010at 
lAharonian & Tayloii2010l) . There are a multitude of potential 
signatures to be explored. 
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